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O ' A dynamical treatment of Markovian diffusion is presented and several applications discussed. 

■ The stochastic interpretation of quantum mechanics is considered within this framework. A model 
OO . for Brownian movement which includes second order quantum effects is derived. 

> 

O '. 1 Introduction 

' Continuous Markov processes, a well developed field of mathematics [Q^IQIj have not yet 
had their full impact on modern physics. The earliest of these, developed by Einstein [0, 
Smoluchowski p], and later refined by Wiener were the original attempts to model diffu- 
sion. They did not have a simple phase space interpretation, and as a result were considered 
Q^l unphysical. The Ornstein-Uhlenbeck model |l^ which did allow a phase space picture, but 
was not a Markov process in coordinate space, became the preferred description. Although 
^ ■ mathematicians pursued the field vigorously, physicists consequently lost interest. The work 
^ . of Nelson |Tl], |T2| was a noteworthy exception. He studied the general problem of diffusion 
and found in a particular model that if a quantity called the mean acceleration were equated 
with the external force, paralleling Newton's equation, then Schrodinger's equation could be 
^ ' derived, suggesting a connection with quantum theory. The results presented here are con- 

■ sistent with this, although a different approach is taken to the same model, and a different 
dynamical assumption is made. 

In section 2 a new technique of analysis is developed for the simple but important case 
of the Wiener process. Velocities as random variables cannot be defined for these processes. 
The sample trajectories of the diffusing particles are nowhere different iable almost surely. A 
Hilbert space approach is presented in which coordinates, velocities, and accelerations are 
noncommuting operators. This is an improvement over an earlier work which introduced 
the concept of ordered expectation [^, and which presented many of the basic ideas of 



the Hilbert space method. The various operators have a physical interpretation as they are 
defined in terms of limits of expectations of random variables. This method is used to derive 
several known results to illustrate consistency. 
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In section 3 the Hilbert space method is apphed to the general Markov diffusion model. 
Operators for coordinates, velocities, and accelerations are defined in the same way as for the 
Wiener process, and they are related to the properties of the underlying stochastic model. 
Commutation rules are derived and equations governing the time evolution of the diffusion 
presented. Section 4 generalizes the technique to three dimensions. The approach is easily 
generalized to any system with a finite number of degrees of freedom. 

The dynamical assumption, introduced in section 5, equates the acceleration operator to 
a sum of external forces plus a stochastic force. This stochastic force is chosen so as not to 
violate momentum and angular momentum conservation in the absence of an external force 
and in the absence of viscosity. The resulting model has four parameters. By varying these, 
models are obtained for Schrodinger's equation and for thermal diffusion. 

It is found that Schrodinger's equation cannot be identified with a unique Markov pro- 
cess, but rather an infinite number of different processes can lead to it. This curious fact 
suggests that the stochastic interpretation of quantum mechanics be modified to incorporate 
this non-uniqueness and perhaps resolve some of the conceptual difficulties which the inter- 
pretation has faced in the past. The non-commutative approach developed here, so similar to 
quantum mechanics, increases the need to explore completely the connection between these 
two theories. 

The thermal model approaches the quantum mechanical Gibbs distribution in the steady 
state limit, with error terms of order fi^, as it should. It can be generalized to include 
viscosity, and it has a Newtonian limit when diffusive fluctuations are small. The usual 
diffusion approximation is found to be a good one under certain circumstances. 

A method for solving the diffusion equations is presented which aids in understanding 
and which utilizes a mathematical similarity between the diffusion theory and the Hamilton- 
Jacobi approach to classical mechanics. 

2 The Wiener process 

A collection of measurable functions w(t), t e (—00,00) on a probability space (Q, E, P) is 
called a Wiener process with zero drift if 



When specialized to the case w(0) = identically it becomes a Gaussian Markov process for 



E{w{t)) = 0, all t. 




r < s,t, 



(1) 
(2) 



t > 0: 



p{w,t) 



1 



-«)2/2l/t 



(3) 

(4) 



V — constant. 




E{^w{tx) X ••• X W{tn)) 

1 1 



^ E{w{t^)w{t2)) X • • • X (5) 



2'»/2 (n/2)! 



TT 
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where tt denotes a permutation of the t subscripts and (^) vanishes for n odd. The Markov 
property is 

P{w{t)eB\w{si),...,w{sn))=P{w{t)eB\w{si)), t = l,n, (6) 

where B is an arbitrary Borel set, and t > Si > Sj+i, all i. 

The density p satisfies the Fokker-Planck or continuity equation 

d 1 

The Markov transition function for this process, expressed as a density, is 

Pt-six, y) = ] exp (-^—yL) , t>s, (8) 



which satisfies the Chapman-Kolmogorov equation 

Psix,y) = J dzPs-uix,z)Puiz,y). (9) 

Joint probability densities may be constructed with the transition functions. For example, 
the joint density for w(t) = y and w{s) = x, t > s is 

p{x,s;y,t) = p{x,s)Pt-six,y). (10) 

It is well known that the sample trajectories for w are nowhere differentiable a.s. As a 
consequence, instantaneous velocity does not exist as a random variable. Therefore define 

^ .■(t + ./2)-u-(i-./2)^ 

e 

which is a random variable for e > 0, but loses this property in the limit e ^ 0. Using eq. 
(^) one finds 

E{w{t)w{t + 5)) = u, 5>e/2, (12) 
E{w{t)w{t -5)) = 0, 5> e/2. (13) 

These equations remain true in the limit e, S provided the inequalities are respected. 
They suggest that a non-commutative approach may be a possible way to study the Wiener 
process. 

Let Ht be the completed real Hilbert space of functions f{x) with inner product 



{f,g) = E{f{w{t))g{wm = J dx p{x,t)f{x)g{x). (14) 
Ht is spanned by the Hermite polynomials. Define a linear operator w on Hf by 

{f,wg)= hm Eif{w{t-6,))wit)g{w{t + 62))), (15) 



where f,g& Hf. This definition is independent of how the Si are taken to zero because the 
expectation in (|15D satisfies 

Eifwg) = J dxdy p{x,t~ 5i;y,t + 52)fix)g{y) 

X E{w{t) \w{t- 6i) = X, w{t + 62)=y). (16) 

This equation follows from the Markov property together with the definition of conditional 
expectation. It is easy to show 

E{w{t) I wit - = X, wit + ^2) =V) = (17) 

Q\ + 02 

and so 

x^ \ I {y — x)^ \ y — X 



(/, wg) = lim J dx dy exp i "^-^^yr ) exp 



2u{t-6i)J MSi + 52)J 5i + 62 

xf{x)g{y)^[{t-6^){t + 62r'/'. (18) 

In this expression t — 61 may be replaced by t (for t > 0) without affecting the limit and the 
resulting expression is a function of 61 + 62 so that the manner in which the limit is taken 
doesn't matter. 

Suppose that g{y) is a polynomial. It has a Taylor's expansion about x 

g{y) = g{x) + {y-x)g'{x) + ---. (19) 
When substituting this expansion into (pISD, only the second term contributes. One finds 

U,^9)={f.^^A- (20) 



dx 



Since / is any element of Ht, if follows that 



wg = ^^^1 9 ^ polynomial. (21) 
If an operator w is defined by simple multiplication, then it follows 

[w, w]g = ug, g a polynomial. (22) 
Similarly, an operator for acceleration may be defined. First define a random variable 



and next define an operator on if, 



(/,%)= lim EU\w{t-5^))w{t)g{w{t + 52))). (24) 



It is easy to show that this vanishes from the fact that 

E{w{t) I w{t - 5i) = X, w{t + 82) =y)=0, 
which follows from eq. (p!?]). Therefore 

w = 0, (25) 

as an operator equation. Likewise all operators for higher time derivatives defined in this 
way are zero. 
Let us define 

H = \w\ (26) 
H is not symmetric when restricted to the polynomials 

{f,Hg) = {Hfp,g/p), (27) 

and it satisfies 

[H,w] = iyw. (28) 

H governs the time evolution of expectations as can be seen by theorems I and II. 
Theorem I For / a polynomial 



d 
dt 



WHt)))= (i,^[^,/H])- (29) 



Proof 



i,-[^,/H]) = (30) 

1 

dx-uf{x)—p{x,t) (31) 

dx/(x)^p(x,t) = ^E(/(«;(t))), (32) 

where the Fokker-Planck equation (|^ has been used in going from (^) to (^). 
The following operators are well defined when acting on polynomials. 

w{t + s) = e(i/^)^"u7 e-(i/'')^" = w + sw. (33) 

They satisfy commutation rules 

[w{ui) ,W{U2)] = l^iUi-U2). (34) 

These operators play a role similar to the coordinate operators in the Heisenberg represen- 
tation of quantum mechanics as can be seen from theorem II. 
Theorem 11 

E{w{ti) X ■ ■ ■ X w{tn)) = {l,w{ti) X ■ ■■ X w{tn)), U < U+i. (35) 
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Proof Let Si = ti — t. Using 



one has 



dx 



Integration by parts of the right hand side, using the Fokker-Planck equation, yields 



d 



d 



dx p{x, t){t + Si)^ — X + S2i'Tr- X ■ ■ ■ X I x + s„z/ 



dx 



dx 



dx 



But, using (HI), this may be written 



(1, wih) X • • • X wiQ) = J2 E{w{h)w{U)) 



1=2 



(36) 



where T, a notational device, denotes time ordering of the product, with later times standing 
to the right. From eq. (|^) straightforward combinatorics yields 



E{w{h) X ■ ■ ■ X witn)) = J2 E{w{h)w{ti))E 



i=2 



\ 



I 



(37) 



Subtracting (|37|) from (0) yields 

n 
1=2 

where 

Dnih, . . . , y = 1^1, T n Witj)j - E{w{ti) X ■ ■ ■ X W{tn)). 

But, Di = D2 = Qhj direct computation, and the theorem follows immediately by transfinite 
induction. 

As a corollary to theorem II, for fi polynomials: 

E{fi{w{ti)) X ■ ■ ■ X fniwitn))) = (1, fl{w{ti)) X ■ ■ ■ X /„(w(t„))), U+i > U 

which follows by expanding the polynomial products using theorem II. 

For the Wiener process, time translations correspond to similarity transformations on 
the coordinates as in (|33D . It is interesting to compare this with quantum mechanics where 
time translations correspond to unitary transformations. This is a fundamental difference. 
For example, a product of different time coordinates in quantum mechanics is not hermitian 
and does not even have a real expectation. 

An ordered expectation may be defined in the following way. Let 



Wi{t) = w{t), W2{t) = W{t), Wsit) = w{t) 
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and define 



Wj^{t) X • • • X w,„(t) = lim E{wj,{ti) X • • • X (38) 

e/(ti+l-ti)^0+ 



Ordered polynomials f(w,w) are therefore defined. Using theorem II it follows 

f{w, w) = (1, f{w,w)) ^ J dx p{x, t)f (^x, . (39) 

Likewise, for ordered polynomials involving w, w, and w one finds 

f{w, w, w) = (1, f{w, w,w)) = J dx p{x, t)f ^x, i/^, j . (40) 

These ordered expectations satisfy a kind of completeness relation. Defining 

f{w,w)5{w-z) = J dxp{x,t)f ^x,iy-^^ S{x - z), (41) 

5(w-z)f(w,w) = j dxp{x,t)5{x-z)f(^:,v^, (42) 

the following results are obtained after several integrations by parts: 

/ / -N f{w, w)5{w - z) d{w - z)g{w, w) = f{w, w)g{w, w). (43) 

This result is similar to quantum mechanics where states of definite coordinate form a com- 
plete set of states for a spinless non-relativistic particle. 

The S function can be used, at least formally, to derive the correct probability densities. 
For example, 

{l,6{w{ti)-y)) = 1 da;p(a;,t)(5(e(i/'^)^(*i-*)a;e-(V'^)^(*i-*) -I/) 
= / dxp(x,i)e(V-W-*)5(x-y)e-(V-W*^-*) 

= J dxS{x-y)e^^^''^''^''-'^p{x,t) 

= P(y,h), H^\i.^^. (44) 



Likewise, 
and generally 



(1, T5{w{U) - y)5{w{t2) - z)) = p{y, h; z, h), (45) 
(1, T5{w, ti)-yi)x---x 5{w(tn) - y„)) = p(yi, h;...; y„, Q. (46) 
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3 Generalized Brownian motion 



The non-commutative approach outhned for the Wiener process has an analog for a deeper 
class of stochastic processes defined formally by the stochastic differential equation: 

dx{t) = b{x{t),t)dt + dw{t). (47) 

It is more precise to consider a stochastic integral equation which is formally equivalent to 
this: 

t 

x{t) - x{0) = I b{x{t'),t')dt' + w{t). (48) 



This equation has been studied extensively in the literature Nelson |11, 12] examined 

the relationship between this theory and non-relativistic quantum mechanics. 

Higher dimensional diffusion will be considered, but first the one dimensional case shall 
be examined. Doob [l|] has shown that ( ^Sf ) is solved by iteration provided b{x, t) satisfies 
certain conditions: it must be a Baire function in the variables x and t, it must satisfy the 
bound b < k\^l + x^, and it must satisfy a uniformity condition: 

\b{xi,t) — b{x2, t)\ < K\x2 — Xi\ 

for some constant K. He also showed that x{t) is a Markov process and that its sample 
trajectories are almost all continuous. 

The Markov transition function plays an important role in the analysis of this process. 
Rather than working with the usual transition function, which is a distribution, it is conve- 
nient to work with its associated density, assuming that this exists. The transition density 
satisfies a Chapman-Kolmogorov equation: 

P{x,t;y,s) = J dz P{x,t : z,u)P{z,u;y, s), t>u>s (49) 

where 

d 

P{x,t;y,s) = —P{x{t) < x\x{s) =y), t > s. (50) 

With regularity assumptions, Doob derives a forward (FE) and backward (BE) equation 
for P: 

FE: ^^p(x,t;y,s) + -^b{x,t)P{x,t;y,s)-^u-^P{x,t;y,s) = 0, (51) 

BE : ^ P(x, t:y,s)+ b{y, t) ^ P(x, t; s) + ^ ^ P(x, t; y, s) = 0, (52) 

where t > s, and u is the diffusion parameter for w in eqs. ( ^Tj ) and (P5|). He also shows 

'x{t + h) - x{t) 



lim E 



h^o+ y h 

and 



x{t) =x \ = b{x,t) (53) 



{x{t + h)- x{t)) 



hm E , 

h-*0+ \ h 



x(t) =x] = u. (54) 
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With these results in mind, let p{x,t) denote the probability density for x(t) at time t 
(assuming that it exists) 

p{x,t) = ^P{x{t)<x). (55) 
Denote by Ht the real Hilbert space of functions on R with inner product 

(/, g) = I dx f{x)gix)pix, t) = E{f{x{t))g{xm- (56) 

R 

It shall be assumed that p{x,t) is a function of fast decrease as a; — »• oo. Therefore, Ht 
contains all of the polynomials, and these are dense in Ht. Define random variables: 

Define an operator x on Ht which has matrix elements for / and g polynomials: 

{f,xg)= hm E{f{x{t-6^)x{t)g{x{t + 62))). (59) 

Using the Markov property, which is the same as (|^), it follows 

E{f{x{t-6i))x{t)g{x{t + 62))) 

dx dy p{x, t - 6i)P{y, t + S2; x,t~ 5i)f{x)g{y) 
xE{x{t) I x{t - 6i) = X, x{t + 82) = y). (60) 



Using the forward and backward equations, the conditional expectation in ([60| ) can be cal- 
culated. The result is that ( |60D can be rewritten as 

= J dxdydz p{x,t - 5i)f{x)g{y)S^P{z,t;x,t - 5i) 

d ^ 
xP{y, t + 62; z, t)b{z, t) + P{z, t;x,t- 6i)u — P{y, t + 62; z,t)j, (61) 

where the e — > limit has already been taken. Since the sample paths are continuous, it 
follows that 

limP(x,t;?/,s) = (5(x - y). (62) 
For the limit 61,62 0, (p2| ) can be substituted into ([6T|). The result is 

{f,xg)=J dxp{x,t)f{x) (^{x,t) + iy-^ g{x) (63) 

which may be rewritten 

{f,xg) = J dxe^-^/(x)z/^^7(a:)e^+^ (64) 
9 



where R and S are defined by 



R 

b 



-ln(p(x,t)), 
.^iR + S) 



(65) 
(66) 



and S is defined only up to an arbitrary additive function of t. 
An operator for acceleration is similarly defined as follows 

(/, xg) = lim E{f{x{t - 6^))x{t)g{x{t + S2))) 

but, 



(67) 



E{f{x{t-6i))x{t)g{x{t + 62))) = J dx dy p{x,t - 6^) 

xP{y, t + 62;x,t- 5i)f{x)g{y)E{i{t) \ x{t - 5i) = x, x{t + ^2) = y). (68) 

Using the forward and backward equations and the definition of conditional expectation, 
this can be rewritten as 

j dxdydz p{x,t - 5i)f{x)g{y)P{z,t; x,t- 5i)P{y,t + §2; z,t) 



db(z,t) 1 d%(z,t) 1 a ,2. 



dt 



Substituting the 6 function of ( |62l) into this expression yields the proper limit for 
finds 

/ f ■■ \ /■ 1 / X N / \ \db(x,t) 1 d%{x,t) 1 d 
{j,xg)= / dx p{x,t)f{x)g{x) ' ■■ 



dt 



—u ■ 
2 



dx"^ 



2 dx^ 



So that 



db 1 d^b Id 



X 



+ 7:^ 



+ 



dt 2 dx^ 2 dx 
or defining a stochastic potential U by 

dU 
dx 



X 



([fl|) becomes 



dU db l^d% \d_^2 
dx dt 2 dx"^ 2 dx 
But using (|66|) this can be transformed into: 



0. 



M+S 



d 



R+S 



To supplement this equation, write the density at time t as 

p{x,t) = J dy p{y,s)P{x,t;y,s), t > s. 



(69) 
One 
(70) 

(71) 

(72) 
(73) 

(74) 
(75) 
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Therefore, p must satisfy the forward equation 



dp ^ d ^ 1 d"^ Q 
dt dx 2 dx"^ 



Using (|65| ) and (|66D this becomes 



dp 
di 



d_ dS_ 

dx dx ^ 



(76) 



(77) 



Equations (^4|) and (iTTf) are easily shown to be equivalent to 

d 



1 2 5' 



(9t 



(7J 



These are the basic form of the diffusion equations which will be used here. They are referred 
to as the Markov wave equations. They may be understood heuristically in the following 
way. Ordered expectations may be defined for x(t), x{t), and x{t) as in eq. (|38|). For / an 
ordered polynomial, one finds 



f{x,x,x) 



1,/ x,x, 



dx 



dxe^ ^ f { X, u—, ——— ] e 



d dW 



R+S 



dx ' dx 



The commutation rules implicit in (p3| ) and ( |64D are 



X, x = z/, 



(79) 



(80) 



and from these one deduces the following rule, using the chain rule of differentiation together 
with the definition of ordered expectations: 



dt 



f{x,x) 



dt 



'lJ{x,x)) = {l,[HJ{x,x)])/v, 



where 



H=-x^ + U. 
2 



But (PTF) is equivalent to 



dxe 



R-S 



d_ 

dx 
1 d^ 



R+S 



-V 



2 (9x2 



+ Uj[x,u 



d_' 

dx 



^R+s 



(81) 
(82) 



(83) 



which must be satisfied for arbitrary /. This is possible only if the Markov wave equations 
(iTSl) are satisfied. 
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The situation is similar to quantum mechanics. The arena of the theory is a Hilbert 
space Ht- Coordinates and velocities are operators on this space with commutation rules 
given by (|80|). There is an operator H which determines time derivatives. 

The Wiener process is a special case of this theory. When 6 = then it is seen from ( ^Sf ) 
that x{t) is a Wiener process, provided a;(0) is identically zero. If h is zero then one may 
chose + 5* = 0. The Markov wave equations become 

which force that f/ = 0, and 

1 3"^ f) 
which is just the diffusion equation (J^). Eq. ( |6^ becomes 

{f,xg) = I dxe'''fix)u—g{x), (86) 

which agrees with eq. (|20|). If 6 = 0, then x = from (|7lD , which is consistent with the 
Wiener result w = [(25)]. Thus the general theory consistently describes the Wiener 
process as it must. 

4 Three dimensions 

In three dimensions, the stochastic differential equation becomes: 

dx = bdt + dw, (87) 

where w^, Wy, and are assumed to be independent and to have the same diffusion constant 
v. 

E{wi{t)wj{s)) = z/(5jj min(t, s). (88) 
The forward and backward equations become vector equations. For t > s one finds 

FE : ^P(a;, t; y, s) + V, ■ b{x, t)P{x, t; y, s) - ^uA,P{x, t; y, s) = 0, (89) 

BE : -^^P{x, t; y, s) + b{y, t) ■ VyP{x, t; s, y) + ^A^P(x, t; y, s) = 0. (90) 

The Hilbert space Ht can be formed and an operator x defined as in the one dimensional 
case. The velocity operator satisfies 

(/, xg) = J d'xp{x, t)f{x) [b {x, t) + uV] g{x). (91) 
Likewise, x can be defined, and one finds: 

X = —b+^Ab+ -Vb"^ + {V X b) X {b + uV) (92) 
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as an operator equation. 

Examining (p2|), one sees that in general x cannot be written as a gradient of a potential 
function U. The acceleration becomes a gradient only if 



V X 6 = 0. 



(93) 



This case shall be considered in this paper. Magnetic forces require the non-vanishing of 
(^), but these shall be deferred to a future work. 
The commutation rules implied by (PT|) are 



[Xi , Xj J 



0, 



[Xij Xj 



] = u6ij, [xi, Xj] = iy{dibj - djk), 



(94) 



and with the condition (|93|) , the last commutator in ( plj) vanishes, so that the different 
velocities commute. 

If (^3]) is satisfied, then one may write 



R 
b 



uV{R + S), 



(95) 
(96) 



where (^B|) is the defining equation for S which is defined only up to an arbitrary additive 
function of t. Eq. (R^) becomes 



x = V 



z/^ {R + S) + K'A {R + S) + (V {R + S)f 



-vu 



(97) 



and this equation also defines the potential U. 
Analogous to eq. ([76D, p must satisfy 



but (|97D and (|9^) are equivalent to the Markov wave equations: 



(98) 



-A + U 



^R±S ^ ^^^gi?±S_ 

ot 



(99) 



Either 



or 



may be used to derive the equation of continuity 

-V ■ vpVS, 



dp 
dt 



(100) 



so that the probability flux may be identified with upVS, and the velocity field of the 
diffusion with uVS. 

In the same sense as in section 3, time derivatives of operators may be calculated 



(101) 
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where 



H = -x'' + U. 
2 



Some useful formulae which may be derived are 

lira E{x{t)) = J d^xp{x,t)iyVS. 
lim{E((if)-^} = 



l,f{x,x)) = J Sxe^-^f{x,uV)e^+^. 



(102) 

(103) 
(104) 
(105) 



The heuristic arguments of section 3 could equally well be applied here, but since they 
are the same, they will not be repeated. The restriction V x 5 = is essentially the same 



as Nelson's V x v = 12 



5 Dynamics 



In this section the Markov wave equations (|99D are considered. They shall be used to model 
a diffusing particle of mass m. Multiplying (RDI) by m yields 



2 



^ -Hits' 



=Fmz/ e 



R±S 



V = mU, 



(106) 



where V has units of energy. Let the diffusing particle be subject to an external potential V. 
Equating V and V does not work. For example if = l/2A;x^, a harmonic oscillator, then 
one expects a steady state solution with p = e^^ and 5* = —At, A a constant. The equations 
become in this case 



1 1 
-mu A H — kx 

2 2 



Xmue 



R 



The solutions to this equation are not normalizable and are unacceptable on physical grounds 
since they do not even vanish at infinity. 



Nelson's dynamical assumption [|Tl|, |T^ equated the mean acceleration to the external 
force, and this led to Schrodinger's equation. His dynamics are a special case of a more 
general dynamical system which shall now be presented. 

Since V and V cannot be equated, it shall be assumed that 



V = V + h, 



(107) 



where h is some function of R and S. Translation, rotation, and Galilean invariance shall 
be assumed in the limit V ^ 0. Symmetry breaking effects such as viscous damping can be 
included later. With these symmetires h must yeild no net force or torque: 



J d^xp{x,t){-'Vh) = J d^xp{x,t) [xidj - Xjdi] h = 0. 



:io8) 



In addition to these constraints, the dynamics of the diffusing particle shall be assumed not 
to depend on the overall normalization of p. That is if p is a solution, then Ap should also 
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be for A = constant. The reason for demanding this is that if p is non vanishing in two 
widely separated regions, but zero elsewhere, then the diffusion in the two regions should be 
independent. This is possible only if the above condition is satisfied. 

Two terms have been found for h which satisfy all these conditions. It is likely that they 
are not unique, but they do lead to an interesting theory. The form for consideration is 



h = cR — Xmu' 



VP ' 



A and c constants. 



(109) 



The reader may verify that all of the conditions are satisfied. 

A viscous damping force can be included by adding a term kS, k, a positive constant, to 
V. The final form for the equations are: 



-mu'^A + V(x) — Xmu^ — 
2 ^ ^ VP 



+ cR + kS 



ot 



These equations have a remarkable property. Defining 



;iio) 



/3 



1 



1-2A' 



111) 



then they are equivalent to 
1 



A mz/^A + cR + kS + V 



1711/ d 



13 dt 



so long as A 7^ |, in which case /3 = cx). 

The case A > | is of particular interest. It yields an imaginary number for (3: 



;ii2) 



A > 



(113) 



Defining a complex "wave function" 



the equations become 



1 . 
--rriTT^ A + cR + kS + V 



2 m 



mu dip 



(114) 



(115) 



together with the complex conjugate of (|115|) . This equation is similar to quantum mechanics. 
In fact, for the special case c = 0, k = 0, with 



mi' /— 

n = -j^ = muv2X - 1, 



then ([TTSD becomes 



1^ 
2 m 



A + V 



ip = ih 



dip 
'dt'- 



;ii6) 



(117) 
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which is Schrodinger's equation. Nelson |12] has already pointed out this connection. His 
dynamical assumption is equivalent to choosing = 1, or A = 1. It is seen from (|116|) 
that this is only one possibility. Any three numbers m, z/, and A which satisfy ( |116D will 
yield Schrodinger's equation. This means there is an infinity of Markov processes, each 
with different z/, which lead to Schrodinger's equation. One can even consider "fictitious" 
processes for which v is imaginary which lead to the same equation, but which do not strictly 
speaking have a physical interpretation other than statistical. It may be that the traditional 
interpretation of quantum mechanics is based on just such a fictitious picture of a real Markov 
process. 

A second application of the case \ > \ concerns the quantum-mechanical Gibbs distri- 
bution. Consider a particle weakly interacting with a thermostat at temperature T. Let 
this particle be subject to an external potential well Vix). The Gibbs distribution for this 
particle is 

oo 
n=l 



(118) 



This expression has an expansion in terms of fi. Landau [|l3 (p. 99) finds the following form 
for R up to fourth order corrections in ^, and up to a normalization constant 



R[x) 



\ In(PG) 



2T 1 24T2 m 



-(VV)^ + 4^-A^l 
— 12i m 



(119) 



Compare this expression with ( |115| ). Set S to zero since the Gibbs distribution is a steady 
state. Consider the following choice of variables: 



Eq. ( |115| ) becomes 



— A + 2TR + V 
6m 



-.R 



0. 



(120) 



(121) 



This equation agrees with (|119|) up to fourth order corrections in h. To see this, rewrite 

(122) 



( |12T|) as 



^-"2T 



V - — (ar + {vRy 

6m ^ 



Now iterate 



Ri 

R2 



V 






2T' 






1 

2f 


V- 


-(^ 

6m ^ 


1 


v + 




2T 


12mT 



(123) 



AV 



24mT2 



(124) 



which agrees with Landau's result (|119|) . This is an encouraging result. Further iteration 
does not lead to the correct fourth order term and so the model is only good for small h. 
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A time dependent model is possible within the framework of the model. Using ( |12(]| ), it 



IS 



1 

A + 2TR + V + kS ij = i 

6m J 

ip = exp ( R + i s\ . 



h dtp 



h 



This equation is more useful if one defines 



S = vS 



(125) 
(126) 

(127) 



so that the flow velocity is VS*. Then (p.25|) becomes 



A + 2TR + V + -S 

6m z/ 



exp i R±i 



m 



n 



s 



,. h d ( .VSm 
±z— ^ — exp R ± I — - — D 



(128) 



The parameter u/k, is called the mobility /i and is easily measured. In steady state flow, 
( p,28| ) becomes 



2 h ^y/P ^ ^rpj^ ^}ls = constant, 



6m y/p 



(V5) 
V • pVS = 0. 
At high T and slow flow, ( |129D becomes 

2TR + V + -S = constant, 
/i 

which yield upon taking a gradient 

Flux = pVS = -pTVp - pS/V. 



(129) 
(130) 

(131) 

(132) 



In most practical applications, eq. ( |132| ) for the probability flux is used together with the 
equation of continuity 



The two together yield 



^ = DAp + pAV, D = pT, 



(133) 



which is the usual diffusion equation used in physics. It is only approximate, and in the 
Markov description it is a kind of quasi-static approximation. When flow velocities are high, 
or when quantum effects are important, then (|133D is invalid. 
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A Newtonian limit is obtained from the Markov equations ( p.28|) when the diffusing par- 
ticle is well localized. The following may be derived by direct computation using (|128|) , 



m^E{xit))= J d=^xp(x,t)(-Vy)--^^E(a;(t)). (134) 



If W changes slowly over the region where p is non vanishing, then if x{t) is the mean 
position of particle, one obtains 

mx !^ —Wix) — —X, (135) 

as a trajectory equation which is of course just Newton's equation. 
In the limit ^ ^ 0, eq. ( p.28|) becomes equivalent to two equations: 

-m(ysy + 2TR + V +-S = -m^, (136) 
^^=-V-pVS. (137) 



To understand these, consider ignoring the second and fourth terms in ( |136| ). It becomes 

(ys) + V = -m§. (138) 

This is essentially the Hamilton- Jacobi equation for classical mechanics. Initial conditions 
on R and S must be specified for a solution. 
As an example, consider the case 

S = V -x, t = 0, V = constant. (139) 

With this initial condition, eq. ( p.38| ) can be solved for 5*. Define x{t, xq) by 

d 

—x{t,Xo)=VS{x{t,xo),t), x{0,Xo)=xo. (140) 
ot 

Then a;(t, Xq) is a solution to Newton's equation 

— mx{t, Xq) = -W{x{t, Xo)), X{0, Xo) = Xq, x{0, Xq) = V. (141) 

The solution to ( |137| ) for p is 

p{x, t)= J d^Xop{xo, 0)S^{x - x{t, Xo)). (142) 

This describes the statistical evolution of an ensemble of Newtonian particles all with initial 
velocity V, but having different starting positions distributed by p{xo,0). 
The viscous term in ( 136 ) can also be included in this manner: 

-mS =-m(vsY + -S + V{x). (143) 
2 ^ p 



Defining Xq) by ( 140 ) one finds 



m 



x{0,Xo) = Xq, x{0,Xo) = V, 



(144) 



and p is again obtained from ( |142|) . The interpretation is the same except that the ensemble 
of Newtonian particles is subject to an additional viscous force. 

The only outstanding term in eq. (|136|) in this Hamilton- Jacobi point of view is the 2TR 
term. It can be included iteratively. Consider again the initial condition S = V 
example. First ignore 2TR, solve ( |134| ) and ( |144| ) , and calculate p from ( |142| ) . The first 
approximation to R is 



Ri 



In 



d^a;op(a;o, 0)6^ {x — x (t, xq)) 



(145) 



This function can next be substituted for R in eq. ( |136| ). Qualitatively this term acts hke 
a repulsive potential forcing the particles away from the center of mass. This causes the 
distribution to spread more than the spreading inherent in the Newtonian ensemble. A 
similar treatment is possible for the general model (|115|) , but his will not be presented here. 

The h ^ limit is equivalent to the limit A = |. As a final point, the case A < | in 
( |112|) leads to unphysical solutions as can be seen by examining the steady state harmonic 
oscillator. 

The model which has been presented offers a possible description of diffusion when quan- 
tum effects cannot be ignored, and to second order accuracy in h. In its present form it 
cannot be made to agree to all orders of h with the Gibbs distribution in the steady state 
limit. Undoubtedly part of the problem is that h in ( |109| ) is too simple a function in the 
present formulation. The only other work in which the author is aware of which attempts 
a description of joint thermal and quantum diffusion is the work of Schwinger 
connection with the present theory is not clear at the present time. 



151, but the 



6 Conclusion 

The non-commutative approach presented here is a natural way to analyze Markov processes 
and diffusion. It dramatically suggests a deep connection between the Markov theory and 
quantum mechanics. The fact that the Markov description of Schrodinger's equation is not 
unique is a new and presumably relevant result for the stochastic interpretation of quantum 
mechanics. The models presented for diffusion with quantum effects, if verified, may have 
practical value. The h —>■ limit of these may, under certain circumstances, be an alternative 
to the Ornstein-Uhlenbeck model of Brownian motion. 
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